Non-Sequential Double Ionization is a Completely Classical Photoelectric Effect 
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We introduce a unified and simplified theory of atomic double ionization. Our results show that at 
high laser intensities (I > 10 14 watts/cm 2 ) purely classical correlation is strong enough to account 
for all of the main features observed in experiments to date. 



Short-pulse lasers with high peak intensities (10 14 < 
I < 10 16 , in watts/cm 2 ) now produce multiphoton gener- 
ation of double ionization, the two-electron photoelectric 
effect, with surprising results. To summarize briefly, the 
experimental data show that two atomic (or molecular) 
outer-shell electrons are highly correlated when photo- 
ejected, with a double ionization rate that can be 1- 
million times higher than uncorrelated sequential theory 
1] allows, so the process is called non-sequential dou- 
ble ionization (NSDI). The first laboratory results were 
reported in 1992 and 1993 [3 0, showing an anoma- 
lously high double ionization yield, the principal ex- 
perimental signature of NSDI. Additional data is be- 
ingreported from momentum spectroscopy experiments 
1111,111111,11113113 The momentum distribution data, 
along with the ion-yield data, serve as the benchmarks 
for various theoretical models. 

The mechanism that makes NSDI correlation so ef- 
fective is far from settled, and theoretical exploration 
has been extensive [13 El Q El El G3 Almost all 
existing calculations refer to, or are closely guided by, 
a single few-step rescattering model 0, Qij, which is 
based on an imagined picture in which one electron es- 
capes the atom by quantum tunneling through a field- 
lowered barrier and is then phase-coherently and clas- 
sically forced by the laser away from and then back to 
the core where a quantum collision liberates both elec- 
trons at once (consistent with the term non-sequential). 
However, the patchwork of ad hoc elements typically em- 
ployed 20] has not been claimed to make a complete, i.e., 
self-contained, theory. It is the purpose of this note to 
show that a self-contained theory exists that is compati- 
ble with essentially all prominent features of NSDI. 

Our theory is dynamically classical, and discards all 
aspects of quantum mechanics including tunneling. It is 
built on the need for strong electron correlation to ex- 
plain NSDI, and so must be intrinsically a two-electron 
theory. We do not advocate such a theory for an electron 
that does not have the advantage of a strongly correlated 
partner. It turns out that entirely classical interactions 
are adequate to generate very strong two-electron corre- 
lation, as observed in NSDI, and quantum theory is not 
needed. Of course atoms are quantum objects but in such 
strong fields as are used for NSDI it is mainly electron 
physics rather than atomic physics that determines the 
experimental outcome. In this sense the early remark of 
Corkum [Tsj advocating the adoption of a plasma per- 
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FIG. 1: Energy vs. time plots showing 4 distinct stages of 
NSDI. The orange and blue lines track the energies of two 
electrons. After repeated exchanges of energy while both are 
bound, one electron (orange) acquires positive energy but re- 
turns a number of times to the core. Its strong recollisions 
show up as spikes in the energy line of the still-bound inner 
electron. After several spikes the phasing is right to liberate 
the inner electron and the doubly ionized pair exhibit free- 
electron jitter motion. 



spective was quite appropriate. 

We note that a new form of energy analysis is very 
helpful. The graph in Fig. shows the sum of kinetic 
energy, electron-nucleus binding energy, e-e correlation 
energy and laser field interaction energy for each of two 
electrons during the laser pulse. Such Newtonian en- 
ergy trajectories universally display the same sequence 
of stages, which we propose to accept as the dynamical 
signature of a high-field double-ionization event. 

For simplicity we have made most of our classical cal- 
culations with a one-dimensional model, but the results 
don't depend strongly on this. The results shown in 
Fig. n were calculated one-dimensionally but a fully 
three-dimensional calculation displays exactly the same 
sequence of stages, as shown in Fig. [21 A quick inspec- 
tion of these energy plots shows that we can refer to the 
stages of high field double ionization as initiation, recol- 
lision, ionization and jitter. This characterization is not 
far from that conjectured by Yudin and Ivanov [TEI ] on the 
basis of selected Newtonian trajectories for a single elec- 
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FIG. 2: Classical high-field energy trajectories calculated in 
three dimensions, showing the same characteristic four stages 
of NSDI, as in Fig. □ 



FIG. 3: The "knee" signature of NSDI data is clearly pre- 
dicted even in a classical calculation. The dashed curve (red) 
is inserted by hand to indicate theprediction made by se- 
quential quantum tunneling theory 



tron. Our pictures show details in a time domain not seen 
before (2e dynamics prior to single ionization), and they 
also reveal a previously unremarked electron-pair phas- 
ing. Below we correlate this phasing with momentum 
properties of the electron-ion products in an intuitively 
appealing way. 

One can easily describe the characteristic features of 
each of the four NSDI stages as follows. 

• Initiation stage: Both electrons are confined in bound 
orbits in the nuclear potential and many rapid e-e inter- 
actions occur. One electron can easily take energy from 
the other, and escape without need for tunneling. 

• Recollision stage: The semi-liberated electron returns 
to the core repeatedly. During each return the efficiency 
of energy transfer rests on the relative motional phase 
between the two electrons. 

• Ionization stage: A final collision leads to highly cor- 
related double ionization, which we find usually occurs 
after several recollisions, not the first. 

• Jitter stage: The two electrons exhibit the jitter oscil- 
lations characteristic of free electrons. The oscillations 
are exactly in phase or out of phase with each other. 

We note that in an entirely classical picture one is 
able to work with a self-contained and almost completely 
unified theory. That is, in contrast to mixed classical- 
quantum patchworks, there are few adjustable parame- 
ters, and ad hoc decisions about timings, cross sections, 
matrix elements, etc., are not needed. Moreover the cal- 
culations are complete in the sense that they begin when 
the field turns on and continue without modification to 
whatever later time is of interest, and are exact in the 
sense that perturbation theory plays no role. A classi- 
cal prediction of double ionization might at first be con- 
sidered accidental, except that Newtonian calculations 
produce the anomalous high-yield "knee" , the principal 
signature of NSDI, as shown in Fig. [3J where the ver- 
tical line shows the location of the intensity threshold 
predicted by the old two-step theory. Moreover, we find 
no strongly correlated double-electron effect that needs 
quantum mechanics in order to be understood, not even 
the initial liberation of the first electron, which in the 
two-step model is assumed to originate in quantum tun- 
neling (remarkably, the need for this Ansatz seems never 



to have been tested for two-electron phenomena) . 

Our classical method is straightforward once the en- 
semble of initial conditions is described (for this, see [2l|). 
We employ a large microcanonical ensemble (100,000 - 
500,000 members) of independent two-electron atoms, 
each with initial total energy equal to the energy of 
the corresponding two-electron quantum ground state. 
Since NSDI is understood as a universal high-field phe- 
nomenon to be found in essentially any multi-electron 
atomic species we try to make non-specific calculations as 
far as possible. One nonspecific component is the familiar 
quasi-Coulomb model for the one-dimensional electron- 
electron and electron- nucleus potentials [22] V(x) = 
— l/\/x 2 + a 2 , or its analog in three-dimensional calcu- 
lations. When we take the soft-core parameter the same, 
a = 1, for all the interactions, the initial two-electron 
energy is fixed at the ground-state energy -2.24 a.u. We 
choose an 8-cycle or 10-cycle (25-30 fs) sinusoidal laser 
pulse with the wavelength 780 nm (frequency u) — 0.0584 
a.u.) and a trapezoidal envelope. By integrating the 
Newtonian equations of motion we can then numerically 
follow any two-electron trajectory and the ensemble of 
them is suited to a statistical analysis that can be com- 
pared with experimental results and with the results of 
other theoretical models. 

The final momentum distribution of NSDI ions of Ne, 
He and Ar atoms measured using COLTRIMS are dif- 
ferent, as shown in panels (a), (b) and (c) in Fig. 21 
The distribution of Ne shows a double-peak structure 
with a valley at zero 0, , the distribution of He shows 
a double-peak structure with its valley nearly filled 
and the distribution of Ar shows a broad peak centered at 
zero O ne can interpret that there are two groups 

of ions or trajectories: one group gives a distribution of 
ions peaked at Zero (Z) momentum, and the other gives 
two distributions of ions peaked at Non-Zero (NZ) mo- 
menta. These two base distributions superpose to give 
the observed distributions. In Ne, the main contribu- 
tion comes from the NZ ions with minimal contribution 
from the Z group ions. In He, the main contribution also 
comes from the NZ ions, but it has more contribution 
from Z ions than in Ne in order to fill the valley between 



3 





-5.0 -2.5 25 5.0 
Ion momentum (a.u.) 



2 4 -4 -2 

Ion Momentum (a.u.) 

FIG. 5: Final ion momentum distributions calculated classi- 
cally using the classical ensemble method. The letters Z and 
NZ label the peaks of the broad but distinct groups of tra- 
jectories with small or Zero ion momentum and substantially 
Non-Zero ion momentum respectively. 



FIG. 4: Final ion momentum distributions of neon helium 
[f| and argon Q measured using COLTRIMS. Panels (b) and 
(d) show the momentum distribution of He ions at different 
laser intensities. Two solid curves (red) and a dashed curve 
(blue) are inserted by hand under the experimental data to 
identify contributions from the two types of NSDI trajectory. 
One type gives distributions peaked at Zero (Z) momentum, 
and the other gives distributions with symmetrically placed 
Non-Zero (NZ) momentum peaks. The relative sizes of the 
Z and NZ peaks in each distribution depend on the atomic 
species and the laser intensities. 

the double peaks. As for Ar, the number of NZ ions is 
significantly less than the Z ions to give the broad peak 
structure centered at zero. This interpretation is con- 
sistent with the analysis by Feuerstein et al. @, who 
separate the distributions of NZ and Z ions in Ar 2+ with 
Z ions as the majority. These two groups of ions were 
previously identified in our classical simulation |2^| . 

A new result obtained from our purely classical treat- 
ment of strongly correlated electron dynamics, as re- 
vealed by the new energy trajectories, is the first in- 
tuitively natural explanation for a main component of 
recoil ion momentum distributions such as are obtained 
experimentally for helium, neon and argon. Figs. n an d 
121 show that the end stage of NSDI finds the two ejected 
electrons exhibiting jitter oscillations that are either ex- 
actly in phase or out of phase. These energy oscillations 
are due primarily to the potential energy from interac- 
tion with the laser field. The in-phasc oscillations have 
the electrons typically escaping the binding potential in 
the same half laser cycle after the final electron-electron 
collision, although some of them can have a lag time of 
an even number of laser cycles. This situation gives rel- 
atively non-zero sum of electron momentum and is de- 
noted above as the non-zero (NZ) recoil case. On the 
other hand, the out of phase events find the second elec- 
tron to be field-ionized (rather than directly collision- 
ionized) in an odd half-cycle after the first one departs. 
The electrons emerge on opposite sides of the nucleus. 
This situation gives zero or relatively small sum of elec- 
tron momentum and is denoted as the zero (Z) recoil 
case. 

For completeness, Fig. [S] shows the momentum dis- 
tributions arising from classical calculations, with peaks 



labelled either Z or NZ. Classical dynamics predicts that 
as the laser intensity decreases the NZ fraction and the 
width of the momentum distribution decrease. These 
two effects have been observed experimentally 0, 0, 
(compare Figs. 0] (b) and (d) to see these effects). In 
addition, our classical calculation finds electrons with en- 
ergy above 2t/ p , where U p = I /4u 2 is the ponderomotive 
energy . The number of these electrons is significantly 
less than the number of electrons with energy below 2U p . 
These two findings are consistent with the experimental 
observations on helium, neon and argon from DiMauro's 
group 0,0]. The details of the analysis will be presented 
elsewhere. 

Before concluding, it is appropriate to explain what is 
gained in a classical picture. The biggest advantage is cal- 
culational, since Newtonian dynamics can be followed in 
full detail, something out of the question for an approach 
via time-dependent two-electron Schrodinger theory. On 
the conceptual side the advantages are also substantial 
and significant. Disconnected ad hoc elements disappear, 
and the theory is unified. We note that in the two-step 
picture and its extensions the need for quantum tunnel- 
ing is an ad hoc assumption that has not been critically 
examined. The justification to patch quantum tunneling 
onto classically forced motion requires another ad hoc 
Ansatz. A third ad hoc element occurs in analyses that 
have a classically returning "outer" electron obey a quan- 
tum collision cross section when it encounters the nucleus 
and an artificially quiescient "inner" electron. A more 
sophisticated example occurs in calculations using a se- 
lection of S-matrix elements guided only by the two-step 
picture. The negative effect here is that the rationale of 
S-matrix theory is abandoned, i.e., the next higher order 
of approximation cannot be identified, let alone calcu- 
lated. The ad hoc elements just mentioned are almost all 
inherited from a one-active-electron approach that has 
worked extremely well in single ionization contexts. The 
work of Lewenstein, et al., [24j clearly demonstrates this. 
However a one-electron theory is not suitable for an in- 
trinsically two-electron phenomenon. 

In summary, we have made an entirely classical study 
of the response of two-electron atoms to intense short- 
pulse laser radiation in the very high-field regime. Our 
purpose is not to find agreement with NSDI data specific 
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to any real atom, but to determine whether it is reason- 
able to assign the strong electron correlation associated 
with NSDI to classical orgins. This is clearly the case 
and the consequences are important because a straight- 
forward and numerically easy route is now opened up for 
undertaking systematic and wide-ranging exploration of 
time-dependent and phase-coherent multi-electron effects 
in a strong radiation field. This is more than academ- 
ically interesting because pictures of electrons in such 
trajectories are being used to guide important initia- 
tives in atomic, molecular and optical physics, includ- 
ing attosecond timing of atomic and molecular processes 
[25|. the generation of controlled intramolecular single- 
electron beam currents [2f| and short-wavelength coher- 
ent radiation [2]j , and the use of strong short-pulse laser 
fields to control electron motion in general 28] . In partic- 
ular, we have shown how central is the role of e-e dynam- 
ical interaction in the ionizing dynamics. The strength of 
the e-e interaction has previously been undervalued, but 
it dynamically facilitates the required energy exchange 



between electrons. Different multiple-recollision channels 
produce trajectories with different final momenta, pro- 
ducing similar or opposite jitter phasings, which underlie 
the experimental NZ and Z momentum distributions, re- 
spectively. Most importantly, we have shown that the 
major signatures of NSDI, at least those that are ob- 
served experimentally, can be understood using classical 
physics, leaving only a minor (or possibly future) role for 
quantum effects to play. It is intriguing to guess that 
experimenters have already been registering some three- 
electron effects with NSDI observations without knowing 
it. 
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